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Abstract 

A three-dimensional, space- frequency model for simulation of inter- 
action in free-electron lasers (FELs) is presented. The model utilizes 
an expansion of the total electromagnetic field (radiation and space- 
charge waves) in terms of transverse eigenmodes of the waveguide, in 
which the field is excited and propagates. The mutual interaction be- 
tween the electron beam and the electromagnetic field is fully described 
by coupled equations, expressing the evolution of mode amplitudes and 
electron beam dynamics. 

Based on the three-dimensional model, a numerical particle simula- 
tion code was developed. A set of coupled-mode excitation equations, 
expressed in the frequency domain, are solved self-consistently with 
the equations of particles motion. Variational numerical methods were 
used to simulate excitation of backward modes. At present, the code 
can simulate free-electron lasers operation in various modes: spon- 
taneous (shot-noise) and self- amplified spontaneous emission (SASE), 
super-radiance and stimulated emission, all in the non-linear Compton 
or Raman regimes. 
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1 Introduction 



Several numerical models have been suggested for three-dimensional simu- 
lation of the FEL operation in the non-linear regime [T]-|lfl|. Unlike a pre- 
viously developed steady-state models, in which the interaction is assumed 
to be at a single frequency (or at discrete frequencies), the approach pre- 
sented in this paper considers a continuum of frequencies, enabling solution 
of non-stationary, wide-band interactions in electron devices operating in the 
linear (small-signal) and non-linear (saturation) regimes. Solution of excita- 
tion equations in the space-frequency domain inherently takes into account 
dispersive effects arising from cavity and beam loading. The model is based 
on a coupled-mode approach expressed in the frequency domain and used 
in the WB3D particle simulation code to calculate the total electromagnetic 
field excited by an electron beam drifting along a waveguide in the presence 
of a wiggler field of FEL. 

The unique features of the present model enable one to solve non-stationary 
interactions taking place in electron devices such as spontaneous and super- 
radiant emissions in a pulsed beam FEL, shown in Fig. [H We employed 
the code to demonstrate a spontaneous and super-radiant emissions excited 
when a bunch of electrons passes through a wiggler of an FEL. Calculations 
of the power and energy spectral distribution in the frequency domain were 
carried out. The temporal field was found by utilizing a procedure of in- 
verse Fourier transformation. Super-radiance in the special limit of 'grazing' 
(where dispersive waveguide effects play a role) was also investigated. 



2 Dynamics of the particles 

The state of the particle i is described by a six-components vector, which 
consists of the particle's position coordinates rj = {xi,yi,Zi) and velocity 
vector Vj. Here {x, y) are the transverse coordinates and z is the longitudinal 
axis of propagation. The velocity of each particle, in the presence of electric 
E(r, t) and magnetic B(r, t) = /iH(r, t) fields, is found from the Lorentz force 
equation: 

^ (7iV.) = - ^ [E(ri, t) + X B(r„ t)] (1) 

where e and m are the electron charge and mass respectively. The fields 
represent the total (DC and AC) forces operating on the particle, and include 
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also the self-field due to space-charge. The Lorentz relativistic factor 7^ of 
each particle is found from the equation for kinetic energy: 

where c is the velocity of light. 

The equations are rewritten, such that the coordinate of the propagation 
axis z becomes the independent variable, by replacing the time derivative 
^ = f^.^. This defines a transformation of variables for each particle, which 
enables one to write the three-dimensional equations of motion in terms of 
z: 

^Ii = _JL±v..E (4) 

dz mc^ Vz^ 

The time it takes a particle to arrive at a position is a function of the 
time toi when the particle entered at 2; = 0, and its instantaneous longitudinal 
velocity Vziiz) along the path of motion: 



3 The driving current 



The distribution of the current in the beam is determined by the position 
and the velocity of the particles: 



(6) 



here qi is the charge of the ith macro particle in the simulation. The Fourier 
transform (in the positive frequency domain) of the current density is given 
by: 

+00 

J(r,/) = 2u{f) J J{r,t)e-^'-f'dt^ 



(7) 
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here u(f) = I ^' j, ^ is the step function. 
1 0, / < 

This Fourier transform of the current ^ is substituted in the following 
excitation equations to find the evolution of the electromagnetic fields. 

4 The electromagnetic field 

The Fourier transform of the transverse component of the total electromag- 
netic field is given at the frequency domain as a superposition of waveguide 
transverse eigenmodes 



Ex(r,/) = Y.{C+,i^J)^~"'"'' + C-,{z,f)e^''^'''}S,Ax,y) 

and the expression for the longitudinal component of the electromagnetic 
field is found to be: 



Where k^^ = \j — {k^^ is the cut-off wave number of mode q) and 

C^g{z,f) and C^g{z, f) are the qth mode's amplitude corresponding to the 
forward and backward waves, respectively. Equations ^ and © describe 
the total transverse and longitudinal electromagnetic field (radiation and 
space-charge waves) [TT] . 

The evolution of the qth mode amplitudes C±q{z,f) is found after sub- 
stitution of the current distribution (0) into the scalar differential excitation 
equation: 



Zq 



J±(r,/) +zJ,(r,/) 



■ £±q{x,y) dx dy 
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(10) 



here 



A/;(/) = 

is the power normalization of mode g, and 



■ z dxdy 



. _ _ j +1, propagating modes 
— ^* I —1, CMif: — off modes 



The total electromagnetic field is found by inverse Fourier transformation 
of dHI) and (H: 



E(r,t) = g?|y E(r,/)e+^-2-^*d/ 
H(r,t) = /H(r,/) e+^'^'f' df 



(11) 



The energy fiux spectral distribution (defined in the positive frequency 
domain / > 0) is given by: 



dWjz) 
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E(r,/)xH*(r,/) -zdxdy 



Propagating 



E [\C^,{zJ)\'-\C^,{zJ)\'\^{U,U)} 

Cut-off 

+ 53{c+,(^,/)C:^(z,/)}53{Ar,(/)}(12) 



5 The Variational Principle 

The solution of the equations Q, ^ and (fTn|) for forward waves is done 
by integrating the equations in the positive z-direction for a given boundary 
conditions at the point z = 0. For backward waves the natural physical 
boundary conditions are given at the end of the interaction region z = L^, 
and the direction of the integration is the negative z-direction. 
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In order to take into account excitation of both forward and backward 
waves, we introduce a variational functional 



dz 



where 



2Ar,(/) 



C+,izJ)a,{zJ) + C.,{zJ)P,iz,f) 



dz 
(13) 



^ J^(r,/) +zJ,(r,/) 



^lAx, y) dx dy 



Jx(r,/) +zJ,(r,/) 



■ ^Iqi^^y) dx dy 



(14) 



SF 



The variational derivative of the above functional is: 







SC^qizJ) ^^^^ + P,{zJ) 



5C+,{zJ) 



dz 

' dC.,{z,f) 

dz 



+ ag{z,f)\ + 



d{C.,{z,f) SC+,{z,f)) 
dz 



dz 
(15) 



For arbitrary variations SC±q{z, f) the functional minimizes (i.e., 6F = 0) 
if and only if the equations (fTfl|) are satisfied, and the boundary term is 



SFb- 
resulting in 



dz 



C^,(0, /) 5C+,(0, /) = C_,(L^, /) 6C+,{L^, f) 



(17) 



This enable solving amplifier scheme in which the boundary conditions 
are C_g(Lu,,/) = and C+q(0,/) = 0, as well as an oscillator configuration 
where the boundary conditions are C_g(0, /) = C_g(L^,/) and C+g(0,/) = 

C+q{Lw^ /)• 
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6 Numerical results 



We shall use the code to investigate super-radiant emission radiated when 
an ultra short e-beam bunch (with duration of 1 pS, much shorter than 
the temporal period of the signal) passes through the wiggler of an FEL 
having operational parameters as given in table [H In this case, the power of 
super-radiant (coherent) emission is much higher than that of the incoherent 
spontaneous emission [T2] . 

Fig.[2lshows two cases of dispersion relations: when the beam energy is set 
to 1.375 MeV, there are two separated intersection points between the beam 
and waveguide dispersion curves, corresponding to the "slow" {vg^ < Vzq) and 
"fast" (vg^ > Vzq) synchronism frequencies 29 GHz and 100 GHz, respectively. 
Lowering the beam energy to 1.066 MeV, results in a single intersection 
(at 44 GHz), where the beam dispersion line is tangential to the waveguide 
dispersion curve {vg = v^q — "grazing limit"). 

The calculated spectral density of energy flux in the case of two well- 
separated solutions is shown in Fig. EK- The spectrum peaks at the two 
synchronism frequencies with main lobe bandwidth of A/12 ~ — - — , where 



'spi,2 ~ 



is the slippage time. The corresponding temporal wave- 

'31,2 

packet (shown in Fig. Elb) consist of two "slow" and "fast" pulses with dura- 
tions equal to the slippage times modulating carriers at their respective syn- 
chronism frequencies. The spectral bandwidth in the case of grazing shown 
in Fig. liK, is determined by dispersive effects of the waveguide taking into 
account by the simulation. The corresponding temporal wavepacket is shown 
in Fig. lib- 
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Table 1: The operational parameters of millimeter wave free-electron maser. 



Accelerator 
Electron beam energy: 
Electron beam current: 
Pulse duration: 

Wiggler 

Magnetic induction: 
Period: 

Number of periods: 
Waveguide 

Rectangular waveguide: 
Mode: 



Ek=1^3 MeV 
/o=l A 
T=l pS 

5^=2000 G 
A^=4.444 cm 

1.01 cm X 0.9005 cm 

TEqi 
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Figure 3: Super-radiant emission from an ultra short bunch: (a) Energy 
spectrum (analytic calculation and numerical simulation are shown by solid 
and dashed lines, respectively); (b) temporal wavepacket. 
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